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Normal form theory is a technique for transforming the ordinary differential equations 

O 

describing nonlinear dynamical systems into certain standard forms. Using a particular 
class of coordinate transformations, one can remove the inessential part of higher-order 
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On 
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nonlinearities. Unlike the closely-related method of averaging, the standard development 
of normal form theory involves several technical assumptions about the allowed classes 



O 

of coordinate transformations (often restricted to homogeneous polynomials). In a recent 



paper [1], the second author considered the equivalence of the methods of averaging and 
of normal forms. The references given there, particularly Chow and Hale [2], should be 
consulted for a full treatment of Lie Transforms. 

In this paper, we relax the restrictions on the transformations allowed. We start with 
the Duffing equation, and show that a singular coordinate transformation can remove the 



^ ■ nonlinearity associated with the usual normal form. We give two interpretations of this 
coordinate transformation, one with a branch cut reminiscent of a Poincare section. We 
then show, when the generating problem is linear and autonomous with diagonal Jordan 
form, that we can remove all nonlinearities order by order using singular coordinate trans- 
formations generated by the solution to the first-order linear partial differential equation 
produced by the Lie Transform method of normal form theory. A companion paper [4] 
discusses these methods in a more general context and treats a specific example with a 
nondiagonal Jordan form for the generating matrix. 
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1. Dufflng's Equation: Removing the Nonlinearities 

The second-order problems modeled by Duffing's equation or van der Pol's equation 
under near-resonance forcing are examples of the behavior that concerns us here. We use 
the notation of the Duffing example in [1] , Eqn. (37) : 

x + uj 2 x + ecx + ehx 3 = sRcos (fit) (1) 

where fl ~ u, e is a small parameter, and the other constants are positive. If we manipulate 
the variables and parameters in a standard way (see[l]) by choosing 

R = fl 2 R, h = fl 2 h, c = fic, uj 2 = fl 2 (1 + ea) , r = fit, y = dx/dr = x' 

and changing to complex state variables 

p 0) = \ + iy) » • • -P = \ - iy) , 

then Eqn. (1) becomes 

f 1 1 h R 1 

p' = -ip + ie l--(a - ic)p - -(a + ic)p - -{p + pf + -(e lT + e~ iT ) j (2) 

and its conjugate. 

In [1] , standard methods are used to show that the coordinate transformation 

p = z + eW (z,z,t) , . . .p = z + eW (z, z, r) (3) 

with 

W = -\{a + ic)z + fe- + \ (V - 3zz 2 - |) (4) 
can transform the p-equation into the normal form 

z > = -iz + eh (z, z,t)+0 (e 2 ) (5) 

where 

f 1 =K 1 z + K 2 z 2 z + K 3 e- t \ 
Ki = ~2^ a ~ lc > ' K2 = — 2~' Ks = T 



It is the z 2 z term that standard normal form theory considers as essential to keep at the 
cubic order for equations of this type. 

Now consider the coordinate transformation from (p, p) to (u, u) given by 

p = u + eV(u, u, t) 

(7) 



u + e 



W (u, u, r) + ^fi(u, u, t) In (|) 



and its conjugate. Explicit computation shows that u (r) satisfies a differential equation 
of the form 

u =-iu + (e 2 ) . (8) 

By using a coordinate transformation with a logarithmic singularity at u = 0, we have 
removed the first-order term in e, including the u 2 u term which normal form theory tells 
us is essential and irremovable. 

A special case of Equation (5) occurs when f 1 = z — z 2 z (note that this does not 
correspond to a Duffing equation with real coefficients). For this special case, the exact 
solution z(t) = r(r) exp (i0(r)) to the nonlinear equation through order e terms can be 
written down for initial conditions r(0) = tq ^ 1, 9(0) = 9q '■ 

r{r) = T l 9(t) = +9 - r. 

>/rg + (l-rg) exp[-2 £ T]' 

Our transformation Equation (7) generates the approximate solution 

r (i)(r) =r + r (l - r 2 )(er), 0™ = 9(t) = 9 - r 

Graphs of both types of functions starting from the same initial conditions z(0) = (1 + i)/2 
are shown in Figures 1 and 2. 
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Duffing to First Order 
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Figure 2 shows the corresponding graph of the approximate solution r^^r). It is clear 
that the approximate solution spirals outward, crossing the limit cycle after a finite time 
To of order 1/e. One can verify that our approximate solution is the first term in the power 
series expansion of the exact solution in the parameter e. Notice that the approximation 
is good to first order in s, but not uniformly accurate in time. 
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Duffing, Approximation 



On a Riemann Sheet 
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Figure 3 shows trajectories for various initial conditions generated by our coordinate 
transformation V. Here we have restricted V to a single Riemann sheet of the logarithm: 
hence the discontinuity at the branch cut along Re[z] < 0, Im[z] = 0. This branch cut is 
what allows the logarithm to "unwrap" the singularity of the Duffing equation. Thus it 
is natural not to continue increasing 9 with time, but to restart the approximate solution 
every time one crosses the negative real axis. This may at first seem strange, as we 
are making a jump in the nonphysical variable u(t). On the other hand, we are forced 
into this to make the coordinate transformation V(u, u, r) a single-valued function: the 
discontinuity in u is needed to make the dynamics of p = u + eV continuous. Our 0(e) 
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solution now systematically approximates a Poincare first- return map T along the negative 
real axis: 

T (1) (r) = r + r(l - r 2 )(2ne) 

giving a systematic approximation for the next intersection T(r) of a curve which crosses 
Im[z] = at r = Re(z) < 0. This interpretation of the dynamics preserves the qualita- 
tive behavior of the original dynamical system, although admittedly does not produce an 
explicit analytical solution for the dynamics. 



2. Solving the Lie Differential Equation 

Why does normal form theory miss this useful transformation? (Alternatively, how 
does normal form theory avoid this nasty singular transformation?) We must look more 
closely at how the transformation of state variables is found by solving the partial differ- 
ential equation produced by the Lie transform process. 

For the case when the state equations have linear generating terms in diagonal Jordan 
form (A = diag(A a ) with no A a = 0, a = 1, 2, . . . , N) , the state equation governing the 
component x a of x will be 

^ = X a x a + efi 1 Ht^) + p^(f^) + ... 

The equation for the corresponding component Wa = W a (t,y) in the transformation 
x = y + eW^ + (e 2 /2)W (2 '> + ... is the Lie equation 

L(W a ) - ^ + E V, ^ - X a W a = /« (t, y). 

For autonomous f^, the d/dt term on the left is dropped from L. In that case, by 
restricting the class of coordinate transformations to homogeneous polynomials, normal 
form theory cannot remove any nonlinearities in /W that lie in the null space of the Lie 
operator . These nonlinear terms therefore comprise the "normal form": all other terms 



are removed by the coordinate transformation. But why should the null space of L (the set 
of functions satisfying L(W) = 0) have anything to do with finding a particular solution 
to the equation L{W) = /? 

To understand this, consider a different linear operator, corresponding to solutions of 
the forced harmonic oscillator: 



There is something special about solutions x(t) for forcing functions in the null space of 
H. If we force at resonance f(t) = sin(t), then / is in the null space of H : H(f) = 0, 
and we see that the particular solution x(t) = t sin(t) is qualitatively different from the 
solution for f(t) = sin(a;£) for other frequencies ui. If we restrict the class of perturbations 
and solutions to finite sums of harmonic waves, then there would be no solution to H(f) = 
sin(t). Normal form theory makes precisely this kind of restriction: by restricting the 
perturbations and solutions to be homogeneous polynomials or in the Hamiltonian case 
to canonical transformations, they have defined away the possibly singular coordinate 
transformations that we study here. 

We now demonstrate, for perturbations of the form fa (t, y) above, that we can find 
solutions W to the Lie operator partial differential equation 



in complete generality by reducing it to an ordinary one. Using the method of character- 
istics, one discovers (see Courant-Hilbert [3], p. 11, for a related transformation) that the 
coordinate transformation from (y, t) to (£,r) given by 



H{x) = x + x = f(t) 




6 = 2/i ; 



^ = (yl /Xf) )/(yl /Xl ), /3 = 2,3,...,iv ; 



r = yiexp(-Ait) 
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will reduce the partial differential equation to the ordinary differential equation 

Ai£i^ = G«(£i, 6, • • • , Cn, r) = y). 

C?l Vol 

where the V a = W a /y a are now to be considered as functions of (£, r) and the G's are equal 
to the functions on the far right evaluated in the new variables. Our solution (7) removing 
the 0(e) term in the Duffing equation, was generated using precisely this method. 

3. Concluding Remarks 

We have studied the Duffing equation using a new approach to the calculation of an 
approximate solution. What about our methods in general? 

Our use of the method of characteristics is perfectly general: the same coordinate 
transformations used in the theory of normal forms can be shown [4] to remove all non- 
linearities if the space of allowed functions is not restricted, fmportant questions remain 
about the nature of the higher order terms in e, and about the estimates of the finite time 
for which the approximate solution is valid. (Here we find results valid to times of order 
1/e, but for the example of [4] times of order f/e^ 1 / 4 ) are found.) 

Our simple interpretation of the resulting logarithmic transformation as a Poincare 
return map gave us a correct qualitative picture of the global dynamics for the Duffing 
equation. We do not have a general formula for analyzing other systems in this way, but 
we find intriguing the implied link between the singularities of the coordinate transfor- 
mations introduced by the method of characteristics and the qualitative structure of the 
corresponding dynamics. 
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